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^3 \ A variational wave function constructed with correlated Hyperspherical Harmonic functions 

^^ . is used to describe the Helium trimer. This system is known to have a deep bound state. In 

t^ ' addition, different potential models predict the existence of a shallow excited state that has been 

identified as an Efimov state. Using the Rayleigh-Ritz variational principle the energies and wave 
functions of both bound states have been obtained by solving a generalized eigenvalue problem. The 
introduction of a suitable correlation factor reduces considerably the dimension of the basis needed 
to accurately describe the structure of the system. The most recent helium-helium interactions have 
been investigated. 
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I. INTRODUCTION 



The system He-He is greatly interesting both from a theoretical and an experimental point of view, and it has 
been object of intense investigations in the last years. Despite its simplicity, it is not easy to determine whether it 
supports or not a bound state. Experimentally, usual spectroscopy techniques are not suitable to its study, and only 
recently [l^H diffraction experiments proved its existence, with a direct measurement of its bond length < R >. Its 
binding energy has been estimated through the relation \Ei,\ « h /4m < R >^ and the s-wave scattering length has 
been estimated as ao ~ 2 < i? >. The most recent values for these quantities have been quoted in Ref. 0] after a 
r^ '. new determination of the bond length by diffraction from a transmission grating. They are < R >= 98 ± 8 a.u., 
CIh] \Eb\ = 1.1 + 0.3/ - 0.2 mK and ao = 197 + 15/ - 34 a.u.. 

'~~', Theoretically difficulties in the description of the ^He2 arise because the He-He interaction results from the sub- 

. ■ traction of the huge energies of the separated atoms, which are only slightly different. Moreover, a 1% decrease of 

^ ' the strength of the interaction makes the system unbound. As a result of the continuous refinement in the past years 

\^ • of both experimental data and electronic structure computational techniques, several potential curves for He-He ap- 

lO l peared in literature. Most of them are presented and compared in an article review by Janzen and Aziz M. However, 

^D due to the vivid interest in ^He2 , and to the difficulties in making a really accurate ab initio potential, newer and 

[_^ more accurate potentials have recently appeared in literature. Among the ones described in Ref. p], the potentials 

called HFDB |^, LM2M2 [Q, TTY [g have been widely used in helium cluster calculations. Furthermore, two more 

up-to-date curves, namely SAPTl and SAPT2 H are now available. The latter is believed to be the most accurate 

characterization of the He-He interaction yet proposed. These potentials are constructed on a completely ab initio 

calculations made by Korona et al. pO| , using infinite order perturbation theory (SAPT) and very large orbital basis 

set. In addition, retarded dipole-dipole dispersion interaction is included over the range — 100000 a.u. (SAPTl), or 

^! the more appropriate 10-100000 a.u.(SAPT2). 

ri ■ All these five interactions support only a single bound state of zero total angular momentum. In Table I we 

^ |, summarize the different characteristics of each potential, as well as salient properties of the associated bound state. 

IL/ ' The latter has been computed solving the two-body Schrodinger equation by means of the Numerov algorithm . We 

.^H ', used the value h'^/m = 43.281307 K a.u.^. ^From the table we can immediately see that the five potentials do not 

j^ ' differ qualitatively among each other, though there is a spreading in the binding energy of the dimer |_Ef,| of 0.51 

H I mK. The SAPT potentials predict the highest binding energies whereas the LM2M2 and TTY predictions are very 

. - - 1 close to each other and are the lowest ones. The differences observed in the binding energy are reflected in the mean 

values of the radius, < R > and < i?^ >, as well as in the scattering length ao. The estimate bond length < R > can 

be directly compared to the experiment and can be also used for an estimation of the binding energy \Eb\ and the 

scattering length ao through the relations given in the first paragraph. Those values are shown in the last two rows of 



Table I and reasonably agree with the estimation of Ref. 



in particular the results obtained with the LM2M2 and 



TTY potential are inside the quoted errors. We can also observe that the system is strongly correlated as its binding 
energy \Ei,\ results from a large cancellation between the kinetic < T > and potential energy < V >. Its spatial 
extension is considerably bigger than the range of the potential, as it is shown in Fig.l where the LM2M2 potential 
and the dimer bound state wave function $d are plotted. Finally, the scattering length ao of the system is bigger 



than the range of the potential by an order of magnitude. All these features characterize the ^He2 as the weakest, as 
well as the biggest, diatomic molecule found in nature so far. Moreover the bound state at practically zero energy 



suggests the possibility of observing an Efimov-like state in the triatomic compound 11 12(|. 

Along with the observations of small clusters of helium atoms, different theoretical methods have been used to 
study the properties of such systems. From the beginning it has been clear that standard techniques could have 
problems to describe those highly correlated structures and, accordingly, more sophisticated methods have been 
applied. In Ref. ||l^ the diffusion Monte Carlo (DMC) method was used to describe the ground state of He molecules 
up to 10 atoms. The ^Hca has been extensively studied by different methods (see Ref. [Q and references therein). 
Theoretically it has been shown that the trimer has a deep ground state of about 126 mK and a single L = excited 
state of about 2 mK. There are not bound states with L > ||l5[. In Refs. |l^Jl3l the L = excited state has been 
studied, in particular looking at those characteristics identifying an Efimov state. In fact this state has the property 
of disappearing when the interaction strength is tuned with a parameter A. For example, the excited state exists in 
the very narrow interval 0.97 ^ A ^ 1.2 for the LM2M2 potential. Though with slightly different values of A, the same 
property holds for the other potentials mentioned above. Therefore the helium trimer gives the unique possibility of 
observing the Efimov effect, as the narrow range in A where the excited state '*IIe3 appears contains the physical case 
A = l. 

In the present paper a set of correlated basis functions is used to describe the ^Hes molecule. The correlated 
hyperspherical harmonic (CHH) basis has been applied successfully in the ground state description of light nuclei [Q . 
Similarly to the cluster of helium, these systems are strongly correlated due to the high repulsion of the nuclcon- 
nucleon potential at short distances. Essentially the method consists in a decomposition of the wave function in 
terms of the hyperspherical harmonic (HH) basis multiplied by a suitable correlation factor which takes into account 
the fact that the probability of finding any pair of atoms at distances smaller than 3 a.u. is practically zero. The 
correlation factor has been taken as product of one dimensional correlation functions f(r) (Jastrow type). In Fig.l 
/(r) is compared to the dimer wave function, showing that both have the same short range behaviour. 

The variational description of the trimer using the CHH basis is twofold. Firstly we would like to evaluate the 
capability of the correlated basis functions to describe a strongly correlated system. Special attention will be given to 



the convergence pattern of the energy for both the ground and excited states. In Ref. 14 1 calculated binding energies 
of the ground and excited states of the trimer obtained by different groups are given in correspondence with different 
interactions. The solutions of the Faddeev equations as well as variational methods and adiabatic approaches have 
been used in those calculations. For the very shallow excited state of the trimer, only few results using the variational 
method have been reported so far, showing the difhculty of describing this state with the required accuracy using 
such a technique. In the present work we will show that it is possible to obtain high precision upper bound estimates 
and wave functions for both the ground and excited states by solving a generalized eigenvalue problem. Moreover, 
a detailed study of the wave function will be performed. In particular the tail of the wave function will be analyzed 
with the extraction of the asymptotic constants. The second motivation of the present work regards the extension of 
the method to describe larger systems. In fact a complete study of the ground state and excited states of the tetramer 
has still to be performed. In this context the variational technique is promising and the present study should be 
considered a first step along this direction. 

The paper is organized as follows. In the next section a discussion of the CHH basis for the systems of three atoms 
is given. The numerical results for the binding energy of the two bound states are given in Sect. 3. Some properties 
of the wave functions and the asymptotic constants are calculated in Sect. 4 whereas the main conclusions as well as 
some perspectives for the extension to larger systems are given in the last section. 

II. CHH BASIS 

In the present study the interaction between three helium atoms is taken as a sum of three pairwise potentials. The 
Hamiltonian of the system will be 

H = T + J2v{i,j) (1) 

where T is the kinetic energy operator and V{i,j) is the He-He interaction that in the present work will be taken as 
one of the potentials mentioned in the previous section. 

Considering the helium atom as a spinless boson, the wave function for three identical spinless bosons can be written 
as a sum of three Faddeev-like amplitudes 

* = V(xi,yi) + ^(x2,y2) + V(x3,y3) (2) 



where the sets of Jacobi coordinates (x^, y^) (i, j,k — 1, 2, 3 cyclic) are: 

Yi = -|(rj +rfc-2r0 
Each i-amplitude has total angular momentum LM and can be decomposed into channels 



(3) 



(4) 



A symmetric wave function requires £a to be even. Moreover £a + La should be even for positive parity states. 
Let us introduced the hyperspherical variables 

Xi=p cos (j)i, yi = psm(f)i 



(5) 



where p is the hyperradius which is symmetric under any permutation of the three particles and (j)i is the hyperangle. 
In terms of the interparticle distances r^j = jr^ — rj| the hyperradius reads: 



p = --^\/ n2 + rh + ni 



(6) 



Using the set of coordinates [p, Hi] = [p, (j)i, Xi, yi], the volume element is dV — p^dpdVti — p^dpsin (pi cos^ (f)id(f)idxidyi 
The two dimensional radial amplitude of eq. (0) is now expanded in terms of the CHH basis 



$a(a;i,yi) 
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where the hyperspherical polynomials are given by |T3] 



(2)pe^.L, 
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(7) 



(8) 



with A/",. 



a normalization factor and P^' a Jacobi polynomial. The quantum number fc is a non negative integer 



related to the grand orbital quantum number K — l^ + L^ + 2k. The product of the hyperspherical polynomial 
defined in eq.dS) times the spherical harmonics coupled to LM in eq.(0) gives a standard three-body hyperspherical 
harmonic (HH) function with defined total angular momentum. 

The other ingredient in the expansion of eq.(|4|) is the correlation factor, taken in the present work of the product 
(Jastrow) type. Its role is to speed the convergence of the expansion describing those configurations in which two 
particles are close to each other. The use of Jastrow correlation factors has a long tradition in the description of 
infinite systems as nuclear matter or liquid helium ||20| as well as in the description of light nuclei [|l8| . The wave 
function describing strongly interacting structures, in which the interaction is highly repulsive at short distances, is 
practically zero when the distance between any pair of particles is smaller than the repulsive core of the potential. 
The correlation factor imposes this behaviour as it can be seen from the specific form of the correlation function /(r) 
given in Fig.l. The short range behaviour of /(r) is governed by the two-body potential whereas its medium and long 
range form is not critical since the structure of the system will be constructed by the HH basis. A simple procedure to 
determine the correlation function for states in which the pair (i,j) is in a relative state with zero angular momentum 
is to solve the following zero-energy Schrodinger-like equation pl|] 



[ — (7rT + -^) 

m or'' r or 



V(r) + W{r)]f{r)^Q, 



(9) 



where V{r) is the He-He interaction used in the Hamiltonian of eq.(|]). The additional term W{r) is included to allow 
the function /(r) to satisfy an appropriate healing condition. It is chosen as 



W{r) = Woexp(— r/7). 



(10) 



The specific value of 7 is not important provided that the ranges of the additional potential W{r) and V{r) are 
comparable |2l|. Hereafter its value has been fixed to 7 = 5 a.u.. The depth Wq is fixed requiring that /(r) -^ 1 for 
values of r greater than the range of the potential V{r). 



The hyperradial functions u"{p) of eq.(R) are taken as a product of a linear combination of Laguerre polynomials 
and an exponential tail: 

«^(p) = ^A^^„L(^)(z)exp(-|) (11) 

ni 

where z = j3p and /3 is a non-linear variational parameter. Let \a^k,m > be a correlated completely symmetric 
element of the expansion basis, where a denotes the angular channels and fc, m are the indices of the hyperspherical 
and Laguerre polynomials, respectively. In terms of the basis elements the wave function (||) results 

*= E AlJa,k,m>. (12) 

The problem is to determine the linear coefficients ^^ ,„. The wave function and energy of the different bound states 
are obtained by solving the following generalized eigenvalue problem 

E Al,„^, < a', k', m'\H - E\a, k, m >= 0. (13) 

a' ,k' ,m' 

The dimension iV of the involved matrices is related to three indices: the number of angular channels N^, the number 
of hyperspherical polynomials per channel K^ and the number of Laguerre polynomials per channel M^. According 
to the Hylleraas-Undheim-MacDonald's theorem |2^j2^, there exists a one-to-one correspondence between the ap- 
proximate energy levels Ei{N) and the exact levels e^ = Ei{oo), the i-th approximate level being an upper bound to 
the i-th exact level. Mathematically the following relations hold: 

E,+i{N + I) > E,{N) > E,{N + 1) (14) 

and 

lim E,{N)^e, (15) 

The implementation of the method in the specific case of the helium trimer in which two bound states are known 
to exist, consists in solving the generalized eigenvalues problem for increasing values of N, until a convergence is 
achieved in the estimates of the ground state Eq and excited state Ei. Moreover, an optimum choice of the non-liner 
parameter /? can be used to improve the pattern of convergence. 

III. BOUND STATE CALCULATIONS 

The generalized eigenvalue problem of eq. (O) can be solved to find bound states of general value of total angular 
momentum LM. Here we are interested in the ground and excited state of the helium trimer both having total 
angular momentum L = 0. In such a case the angular dependence of each i-amplitude of the wave function reduces to 
a Legendre polynomial Pi{fii) with fii = Xi ■ iji. Moreover, the angular channel with £a = La = is, by far, the most 
important and it has been the first one to be considered. Contributions from successive channels, with ia = L^ > 0, 
arc highly suppressed due to centrifugal barrier considerations and can be safely disregarded as it will be discussed 
latter. 

The matrix elements defined in eq. (O) have been obtained numerically. In general, as the dimension of the matrices 
increases, numerical problems could arise from integrals containing polynomials of high degree. In fact, a high number 
of basis functions is expected in order to describe simultaneously both the ground and excited state, which have a 
completely different spatial extent. On the other hand, the correlation functions introduce a complicated structure 
at short distances. Therefore, a dense grid of integration points is necessary. The integrals have been performed in 
the set of coordinates [p, (/)3,/i3] using a Gauss formula in the variable /ia and a Chebyshev Lobatto formula in the 
variable cos(203). Grids of 300 points for the first case and 3000 for the second have been used. In the variable p the 
integrals have been performed on a scaled grid: 

P" = '' , ^ . (16) 



with the choice h = .07 a.u., x = 1-008, and n„iax ~ 800, covering the range — 5000 a.u.. A numerical accuracy of 
10"'^ niK has been obtained in the calculation of the binding energies. 

The convergence of the eigenvalues has been studied increasing the number of basis elements, restricting the dis- 
cussion to one channel, namely the ia = La = channel. In this case a totally symmetric wave function can be 
constructed for values of the quantum number fc = 0, 2, ...., kmax (no symmetric function exists for k — 1). Therefore, 
the number of hyperspherical polynomials Kq included in a specific calculation is kmax, except for k„iax ~ which 
corresponds Kq — 1. The number of Laguerre polynomials is Mq — mmax + 1 with rrimax the maximum degree 
considered. The total dimension of the problem to be solved is A^ = Kq ■ Mq. 

In Table II the convergence of Eq and Ei is shown as a function of kmax^'m-max for the LM2M2 potential. We 
observe that, while the ground state energy Eq converges with a rather small basis set, for the excited state Ei 
it is needed a much bigger basis (about one order of magnitude bigger). The ground state converged to the value 
_Eo = —126.36 mK with kmax — 20 and rumax = 20 whereas the excited state converged to Ei = —2.27 mK with 
kmax — 80 and rrimax = 32. In order to speed the convergence with respect to the Laguerre polynomials the value of 
the non linear parameter /3 has been optimized. For the ground and excited state we have used (3 — 0.40 a.u.^^ and 
P = 0.10 a.u.~^ respectively. 

After the convergence of the first channel is achieved, the contribution of the channel £a = La = 2 can be evaluated. 
The first four linearly independent totally symmetric basis elements belonging to the second channel correspond to 
values of the grand angular quantum number K = 12, 16, 18, 20. The inclusion of these elements gives extremely small 
contributions to the binding energy of the two states and does not change the estimates given above. It is important 
to notice that the Jastrow correlation factor introduced contributions from channels with ia ~ La > 2 already in a 
calculation limited to the first channel. 

Special attention has been given to the study of the convergence with the non-linear parameter /3. In Fig. 2 we 
reported the ground and excited state energy curves as a function of /3 for increasing values of kmax- The number of 
Laguerre polynomials has been kept fixed at rrimax = 24. For the ground state energy the upper curve corresponds 
to kmax — 0, i.e. only one hyperspherical polynomial has been taken into account, and the lower curve corresponds 
to kmax — 20. Results with larger values of kmax are not shown since they completely overlap with the result at 
kmax — 20. For the excited state energy the different curves correspond to kmax = 20,40,60,80. We observed that 
there is a region where the variation of f3 does not appreciably affect the binding energies. 

The variational method provides, in addition to an upper bound to the exact energy of the states, a variational 
estimate of the corresponding wave functions. Through the wave function it is possible to calculate a certain number 
of mean values characterizing the ground and excited state. In Tables IITIV, we computed the binding energy, the 
mean value of the kinetic energy, the potential energy, the interparticle distance rij , and the distance Vi between the 
j-particle and the center of mass. The HFDB, LM2M2, TTY, SAPTl and SAPT2 interactions have been considered. 
Other than typical observables, we also computed the asymptotic normalization constants Ci, for the two bound states, 
as defined in Refs. p4,M and briefly described below, and the percentage P^^ of dimer-like structures in the trimer 
wave functions ^^ (i' = 0, 1). 

The variational wave function is constructed in the present work as a sum of correlated products of polynomials with 
an exponential tail in the hyperradius. The short range behaviour of the wave function is governed by the correlation 
factor whereas in the medium and asymptotic region the expansion should reproduce the structure of the system. It 
is interesting to evaluate the flexibility of the correlated basis to reproduce correctly the asymptotic behaviour. Let 
us introduce the asymptotic function Q^{y): 

e,(y) = 72^^^ (17) 



where g^ = y m \E^ — Ed\/fi , and E^^ and Ed are respectively the trimer (i/ = 0, 1) and the dimer binding energies. 
Here y denotes the second Jacobi coordinate defined in eq.(S). In the configuration in which one atom is far from the 
other two the trimer w.f. for the ground and excited state '^i, behaves asymptotically like: 

*, "^^ c,e,(y)roo(y)$<i(x) (18) 

where $rf is the dimer w.f., and c^ the asymptotic constant. Therefore, the overlap function 

0,{y) = f ^*d{^)Y*o{y)'i'.dyd^ (19) 



is proportional to e *"^/y as y ^ oo. In Fig. 3 we plot the overlap functions 0,^{y) and the asymptotic functions 
Qv{y) for both the ground and excited state. From the figure it is clear that the two curves approach to each other as 



the distance y increases. We also observe the very large extension of the excited state. The asymptotic constants c^ 
are obtained evaluating the ratio Ov{y)/Q^{y) at large y values and are given in Tables III and IV for the different 
interactions considered. For the ground state the five interactions produce similar values of cq, though the values of 
LM2M2 and TTY are slightly smaller. Conversely, for the excited state the result obtained with SAPT2 is smaller 
than that obtained with the other potentials. 

The percentage Pd of a dimcr-like structure inside the trimer is defined as: 

/•OO 

Pd= PMfy'dy. (20) 



The results for the ground and excited state are collected in Tables III and IV for the different interactions. For 
the trimer ground state the probability of a dimer-like structure is around 0.35, whereas for the excited state this 
probability increases up to 0.75 for LM2L2 and TTY interactions. The two SAPT interactions predict a lower dimer- 
like structure, in particular SAPT2. This behaviour is related to the slightly tighter binding predicted by the SAPT 
potentials for the two bound states. As a general remark, the very high value of P\ (nearly 70%, compared to 35% 
of the ground state) suggests that the excited state of ''Hes can effectively be pictured as a third particle orbiting 
around a two particle structure. 

In Fig. 4 we plotted some distribution functions relevant to understand the structure of the two bound states of the 
trimer. Namely, we plotted the pair distribution function p(rij), which represents the probability to find the particles 
i and j at distance r^ , and the mass distribution function m{ri), which is related to the probability to find the particle 
i at distance rt from the center of mass of the system. 

Our results for the ground state agree quite well with the ones published in literature ||l^,|l^,^,^,g7| . For the SAPT 
potentials we find that the bond becomes slightly tighter (« 5%), as a result of the its more attractive well. But 
basically there are not qualitative differences in describing the trimer with any of the different pairwise interactions. 
There is a discrepancy in literature whether the main spatial arrangement of the three particles in the ground state is 
either a quasilinear or equilateral configuration P9[ . Our results seem to agree with the latter. In fact, we can try to 
discriminate between the two by looking at the pair distribution and the mass distribution functions. As it is shown 
in Fig. 4, the probability to find any particle in proximity of the CM is almost zero. This strongly suggests that the 
most probable configuration is the equilateral one. Moreover, the ratio between < rij > / < r^ > is very close to 
the ideal V3 of equilateral triangle for all the potentials we used. Regarding the excited state, the presence of a two 
peaks structure in the pair distribution function agrees with the interpretation of such a state as composed by a two 
particle core surrounded by the third atom at a larger distance. 

For the excited state the results do not depend qualitatively on the potential we use. This may look very surprising 
at a first glance, because this state is suspected to be an Efimov state, and consequently it is expected to be strongly 
affected by any minimum variation of the pairwise interaction. Following Refs. |lq,0,E7|, we studied the behaviour 
of ''He3 as a function of the strength of the pairwise interaction. In Fig. 5 we plotted the energy difference Ei — Ed of 
the system as a function of the parameter A defined by: 

VHe-He - XVl (21) 

where Vl is the LM2M2 potential. We found that this state disappears both increasing and decreasing A, in agreement 
with the claim that it is an Efimov state. Our results agree quite well with Refs. |19,0, where the peculiarity of such 
a behaviour has been widely discussed. To summarize it, we observed that the trimer begins to support an excited 
state at A « 0.975; then, increasing A, the binding energy firstly increases, until it achieves its maximum at A ~ 1.05, 
and successively decreases, until it dissociates. 

In order to compare the different pairwise interactions we assume that due to the very large extension of the w.f. 
compared to the range of the potential, the particles are not sensitive to the particular shape of it, but somehow to 
its average strength. Accordingly we define 



Oi 
C^^' 



(22) 



with x^ L,T, H, SI, S2 in accordance with the LM2M2, TTY, HFDB, SAPTl and SAPT2 interactions and a^ is the 
interparticle distance where the considered potential changes sign, i.e. Vx{(Tx) — 0. In the smaller frame of Fig. 5 we 
reported the different values of A^^ . It is worth to observe that the potentials do not differ so much to show dissimilar 
results, as all the points lie in a small interval of Al = 1. In fact a plot in function of A for the other interactions 
shows that in all cases the physical case A = 1 is on the left of the minimum of the curve, as for the LM2M2 potential. 
The most peculiar feature of an Efimov like system is that it disappears tightening the interaction among its 
components. Physically, such a behaviour could be explained by picturing the system like composed by a third 



particle orbiting around a two particle sub-system. Increasing the strength of the pairwise interaction makes the two 
particles tighter to each other, and the third one evaporates as a result of its very weak bound. 

IV. CONCLUSIONS 

In the present work the helium trimer has been investigated using the most recent helium-helium potential models. 
The helium trimer wave function has been expanded in terms of the CHH basis. Then, the energies and wave 
functions of the ground and excited state have been obtained by solving a generalized eigenvalue problem. The 
Hylleraas-Undheim-MacDonald's theorem assures that the obtained results for the energy of the levels represent 
upper bounds to the exact values. 

The strong repulsion of the He-He potential at short distances engenders some difficulties in the description of 
the three atoms system in terms of an expansion basis. Very large bases are then necessary in order to obtain a 
satisfactory description of the structure of both bound states. The structure is such that the probability of finding 
two atoms at short distances is close to zero and this type of behaviour is difficult to describe using, for example, 
a polynomial expansion. Correlation factors naturally introduce this behaviour accelerating the rate of convergence 
of the expansion basis. In particular the CHH basis combines a Jastrow correlation factor with the HH basis. The 
CHH basis has been used before in the description of nuclear systems |18] in order to take into account the strong 
repulsion of the nucleon-nucleon interaction at short distances. Here, the CHH basis has been used to study five 
different interactions in the description of the trimer. The pattern of convergence for the bound and excited state 
has been studied by increasing the number of basis elements. With a sufRcient number of elements, the dependence 
on the non linear parameter (3 is smooth. Therefore it is possible to obtain a simultaneous description of the bound 
states with high accuracy. The results are collected in Tables III-IV and are in close agreement with previous results 
obtained by different groups using the HFDB, LM2M2 and TTY interactions. The estimates for the binding energy 
are upper bounds to the exact levels and show that the variational method can be used to describe strongly correlated 
systems, as helium trimer, with results that are believed to be among the most accurate ones at present. 

Some interesting aspects of the wave function have been studied. Its asymptotic behaviour in a configuration where 
one atom is moving away from the other two is given in eq.([l8|). In Fig. 3 this behaviour is shown for the two bound 
states. From this study the asymptotic constants Cu have been extracted. In some particular systems the asymptotic 
constants can be measured [pO| . Moreover, the probability Pd of a dimer-like structure inside the trimer has been 
calculated. This quantity gives a clear idea of the spatial structure of the molecule. For the ground state we obtained 
P° « 0.35 whereas for the excited state Pj w 0.70. This latter result suggests a configuration of two atoms in a 
dimer-like state with a third atom orbiting. 

The present work should be consider as a first step in the use of the variational technique with correlated basis 
functions for describing small helium clusters. The extension of the method to study larger systems is feasible. The 
study of the bound states of the tetramer is at present underway and will be the subject of a forthcoming paper. 
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HFDB 


LM2M2 


TTY 


SAPTl 


SAPT2 


Ref. [| 


e 


(K) 


10.95 


10.97 


10.98 


11.05 


11.06 




Vm 


(a.u.) 


5.599 


5.611 


5.616 


5.603 


5.602 




a 


(a.u.) 


4.983 


4.992 


5.000 


4.987 


4.987 




\Et\ 


(mK) 


1.685 


1.303 


1.313 


1.733 


1.816 




(T> 


(mK) 


112.2 


99.43 


99.93 


115.0 


117.8 




~{V) 


(mK) 


113.9 


100.7 


101.2 


116.8 


119.6 




{R} 


(a.u.) 


87.81 


97.96 


97.62 


86.04 


84.24 


98 ±8 


V(R') 


(a.u.) 


119.0 


132.9 


132.5 


116.5 


114.0 




a 


(a.u.) 


170.5 


191.4 


190.7 


166.9 


163.2 




R74m(i?)^ 


(mK) 


1.403 


1.127 


1.136 


1.462 


1.525 


1.1-^0.3/ -0.2 


2{R) 


(a.u.) 


175.6 


195.9 


195.2 


172.1 


168.5 


197 + 15/ - 34 



TABLE I. Characteristic values of the different potentials and their relative bound states. R represents the He-He distance, 
e is the strength of the potential at its point of minimum rm and a is the distance at which the potential changes sign. 





ground state | 


i^max 


5 


10 


15 


20 1 


TTimax 


Eo{mK) 


-Eo(mK) 


Eo{mK) 


_Eo(mK) 


To(mK) 


4 


-120.376 


-120.689 


-120.726 


-120.737 


1769.481 


8 


-126.080 


-126.274 


-126.286 


-126.288 


1662.829 


12 


-126.143 


-126.337 


-126.348 


-126.349 


1660.232 


16 


-126.148 


-126.342 


-126.353 


-126.354 


1660.185 


20 


-126.149 


-126.343 


-126.354 


-126.355 


1660.186 


24 


-126.149 


-126.343 


-126.354 


-126.355 


1660.187 




excited state | 


f^max 


20 


40 


60 


80 1 


TT^max 


Ei{mK) 


Si(mK) 


Si(mK) 


£i(mK) 


ri(mK) 


8 


-1.168 


-1.579 


-1.612 


-1.622 


139.535 


12 


-1.523 


-2.097 


-2.148 


-2.160 


127.139 


16 


-1.555 


-2.150 


-2.222 


-2.237 


123.526 


20 


-1.562 


-2.157 


-2.237 


-2.257 


122.247 


24 


-1.565 


-2.160 


-2.240 


-2.262 


121.943 


28 


-1.567 


-2.161 


-2.241 


-2.264 


121.927 


32 


-1.567 


-2.162 


-2.242 


-2.265 


121.935 



TABLE n. Convergence of the LM2M2 ground state energy Eo and excited state energy Ei for increasing values of the order 
of hyperspherical polynomials k,nax and Laguerre polynomials rrimax- In the last column the convergence for the kinetic energy 
is shown. Basis states with £a = La = have been considered. The non linear parameter (3 has been fixed to 0.40 a.u."'^ for 
the ground state and 0.10 a.u.~^ for the excited state. 







HFDB 


LM2M2 


TTY 


SAPTl 


SAPT2 
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(mK) 


133.0 


126.4 


126.4 


133.8 


135.1 


<T> 


(mK) 


1698 


1660 


1662 


1707 


1715 


- <V > 


(mK) 


1831 


1787 


1788 


1841 


1850 


<ri> 


(a.u.) 


10.38 


10.51 


10.49 


10.36 


10.25 


V<rf> 


(a.u.) 


12.11 


12.28 


12.26 


12.09 


12.03 


< rij > 


(a.u.) 


17.98 


18.21 


18.16 


17.95 


17.77 


V<'-tj> 


(a.u.) 


20.97 


20.71 


20.71 


20.95 


20.84 


Co 




1.22 


1.17 


1.18 


1.24 


1.25 


PS 




0.3614 


0.3310 


0.3311 


0.3619 


0.3539 



TABLE III. Binding energy and mean values of the kinetic and potential energy of the helium trimer ground state calculated 
for different pairwise interactions. The mean values and square root mean values of the distance rt of particle i from the CM, 
and the interparticle distance rtj are also given. In the last two rows the asymptotic constant cq and the probability of a 
dimer-like structure are reported. 
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LM2M2 
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SAPTl 


SAPT2 
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(mK) 


2.735 


2.265 


2.277 


2.788 


2.885 


<T > 


(mK) 


134.1 


121.9 


122.4 


135.7 


137.8 


- <V > 


(mK) 


136.8 


124.2 


124.7 


138.5 


141.7 


<ri> 


(a.u.) 
(a.u.) 


87.24 
103.5 


94.00 
111.1 


93.67 
109.5 


83.07 
96.67 


76.22 


V< ri > 


81.00 


<n, > 


(a.u.) 
(a.u.) 


145.6 
177.7 


157.0 
192.5 


150.7 
183.7 


139.1 
167.5 


125.0 


V< '■I > 


141.1 


Cl 




1.20 


1.24 


1.24 


1.19 


1.09 


p^ 




0.7266 


0.7462 


0.7461 


0.7030 


0.5816 



TABLE IV. Binding energy and mean values of the kinetic and potential energy of the helium trimer excited state calculated 
for different pairwise interactions. The mean values and square root mean values of the distance r^ of particle i from the CM, 
and the interparticle distance rtj are also given. In the last two rows the asymptotic constant cq and the probability of a 
dimer-like structure are reported. 



Figure Captions 

Figure 1. The LM2M2 potential V, its corresponding ground state wave function ^d as well as the correlation function 

/ as a function of the interparticle distance r. 

Figure 2. The ground state energy Eq and excited state energy Ei as a function of the non hnear parameter (3. For 

the ground state the curves corresponding to kmax = 0, 20 are shown, whereas for the excited state the curves are 

given in correspondence to kmax = 20, 40, 60, 80. The number of Laguerre polynomials has been fixed to rrimax ~ 24 

in all cases. 

Figure 3. Overlap functions Oij{y) (solid line) and asymptotic functions Oi/(y) (dashed line) for the ground state 

{v = 0) and excited state (i^ = 1) . 

Figure 4. Distribution functions p{rij) and m{ri) for the ground and excited state of the helium trimer. 

Figure 5. The energy difference Ei — Ed as a function of A. In the small frame the positions in the curve of the 

different values for A^; calculated as explained in the text are given. 
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